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ABSTRACT: This paper is aimed at analyzing the effects of seismic load on piles supporting the electric 
transmission towers erected between Ethiopia and Kenya. Due to the significance of the project, 
failure of the foundations shall be avoided and analyzing the lateral response of piles is crucial. To 
carry out the dynamic analyses, ground motion analyses were executed and a peak ground 
acceleration of 0.13 g has been obtained from seismic hazard analysis for design basis earthquake. 
Acceleration-time history data has been obtained from Pacific Earthquake Engineering Research 
Center ground motion database. Three dimensional finite element analyses were performed to study 
the response of a single pile under earthquake loading by incorporating infinite element boundary 
condition to simulate damping. The non-linear Finite Element Analyses performed by modeling the 
soil using elasto-plastic material model of the Mohr-Coulomb failure criteria indicated that the 
calculated maximum deformations of the unreinforced pile section are not within recommended limits 
of the literature, and the use of 0.6 m diameter pile calls for further detailed analyses. 
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INTRODUCTION 


Transmission tower is an important part of 
power transmission, which provides a safe and 
reliable operation of the electric power 
network. The failure of a power transmission 
tower under seismic actions will affect the 
supply of electricity, causing huge economic 
losses and secondary disasters. When making 
the design of transmission tower, the major 
loads considered are static load, wind load, bolt 
load; while the effects of seismic load are rarely 
considered [1]. Collapse of transmission towers 
is common during an earthquake. In the United 
States, the Landers earthquake in 1992 and the 
Northridge earthquake severely damaged the 
transmission system [2]. The 1955 Kobe 
earthquake in Japan resulted in the destruction 
of large number of transmission towers where 
the major failures were subsidence of the 
foundation and tilting of the tower [3]. The 
1999 Chi-Chi earthquake in Taiwan resulted in 
collapse and distortion on of a large area of a 
tower structure [4]. During the 2008 Wenchuan 
earthquake, more than 20 towers of the 110kV 
transmission line collapsed. Surveys show that 
most damages of transmission towers occur at 


the foundation of the structure [5]. 

Ethiopia is susceptible to two types of seismic 
hazards; namely, earthquake and _ volcanic 
eruptions. [6] cited EM-DAT database that there 
were a total of seven major earthquakes 
between 1900 and 2013, killing a total of 24 
people, affecting 585 people and causing more 
than 7 million US$ economic loses. The 
significant ones are the 1961 Kara Kore 
earthquake (M=7), the 1969 Sardo earthquake 
(M=6.4), the 1983 Wendo Genet earthquake, the 
1985 Langano earthquake and the 1989 
Dobigraben earthquake (M=6.2); some of which 
were fatal [7]. 

Most earthquakes affect tens of square 
kilometers while the severest can destroy areas 
of 2500 km2 or more [7]. The converter station 
is located near the edge of the western plateau 
of the Main Ethiopian Rift (MER), very close to 
the foot of Mt. Damota from the western side. 
Thus, performing dynamic analysis for the pile 
foundation of the electric transmission tower is 
very important. 

The response of a pile subjected to seismic 
excitation depends on the pile itself, mechanical 
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properties of the surrounding soil, the 
interaction between superstructure, pile, and 
soil and the seismic input. The foundation pile 
under investigation is exposed to kinematic 
loading and inertial loading due to earthquake. 
The deformation of the pile is assumed to be 
prevailing because of the stiffness difference 
between the soil and the pile. Secondly, inertial 
loadings are induced in the piles because of the 
accelerations generated within the structure by 
earthquake. 

Consideration is generally confined to lateral 
inertial forces and moments, which are assumed 
to be applied at the pile heads [8]. Novak[9] 
used the generalized Winkler model to simulate 
the dynamic soil-pile interaction system by 
employing linear elastic interaction between the 
two materials. The analytical approach is used 
establish the dimensionless parameters of the 
problem and obtain closed-form formulas for 
pile stiffness and damping. Bentley & Naggar 
[10] studied the effects of kinematic interaction 
on the input motion at the foundation level. 
They included pile-soil separation, slippage, soil 
plasticity, and 3-D wave propagation. [11] used 
a 3D finite element model to obtain the pile 
response under seismic excitation and a load 
applied to the pile cap and column base, with 
consideration of the nonlinearity effects of the 
soil and pile-soil interaction in the dynamic 
behavior of a single pile and pile groups. 
Dongmei &Truman, [12] used a 3D finite element 
model to compute seismic response of a soil-pile 
foundation system, nonlinear analyses are 
performed in the time domain and an infinite 
element boundary condition is used to simulate 
radiation damping. Liuet al., [5] used finite 
element method to compute the response of 
pile-soil-transmission tower under seismic load 
by employing time history analysis. Based on 
the experience from the previous studies, the 
current research has also employed the finite 
element software ABAQUS for the full dynamic 
analyses after determining the seismic hazard 
analyses. The research has thus focused on the 
investigation of the effects of seismic load on 
piles supporting the electric transmission 
towers between Ethiopia and Kenya. 


2. MATERIALS AND METHODS 
2.1. Seismic Hazard Analysis 


The electric transmission tower converter 
station under consideration is located between 
6.9°-7° E and 37.72°-38° N, at the edge the 
western plateau of the Main Ethiopian Rift 
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(MER) that is susceptible to seismic actions. To 
carry out the dynamic analysis, acceleration 
time history record is required although there is 
no recorded data in Ethiopia [13]. Probabilistic 
seismic hazard analysis and deaggregation have 
been performed for the dynamic analysis by 
using time history from previous earthquake 
database. 


2.1.1. Probabilistic Seismic Hazard Analysis 


Earthquake size directly related to the amount 
of energy released in them. The Gutenberg- 
Richter magnitude recurrence model better 
represents earthquake magnitude recurrence 
[14] is given by Equation 1. 

log(/,,) =a-—bM (1) 


m 


where /,, is the rate of earthquakes with 


magnitudes greater than M. The constant a and 
b values are estimated from historical 
observations. Ayele [15] gave a constant b value 
of 0.77, earthquake magnitude value of 7.2 and 
the rate of earthquake 0.593. 


In Chiou and Youngs [19] attenuation models, 
the effect of depth of faulting included the 
depth-to-top of rupture and down-dOip rupture 
width shown in Figure 1. This depth is 
measured to the top of the specific earthquake 
rupture by the empirical equation of [16]. The 
down-dip rupture width is determined by using 
the logarithmic relationships [17]. Style of 
faulting is a useful parameter to include in 
ground motion prediction equation. There are 
three general types of styles of faulting; namely 
reverse, strike-slip, and normal types. The PGA 
calculation has been done using normal slip 
faults because it is the dominant fault 
mechanism [18]. 


Surface 
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Figure 1: Illustration of earthquake source and 
distance measures using a vertical cross 
section through a fault rupture plane. 
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In Chiou and Young attenuation models, 
several distance definition have been used such 

(Rrup), the shortest horizontal distance to the 
surface of projection of the rupture area (also 
referred to as Joyner-Boore distance),Rjg; the 
horizontal distance to surface projection of the 
top of rupture plane, measured perpendicular to 
the fault strike, Ry. Both Ryg and Rx are used to 
model the attenuation of hanging wall effects. 
The distance to the fault is defined to be reliable 
with the specific attenuation relationship used 
to calculate the ground motion. 


The amplitude of the ground motion depends 
on the properties of the material through which 
waves propagate. Site classifications are based 
on the shear wave velocity in the top 30 m 
primary characteristic in the NGA models. Chiou 
and Young [15] add sediment thickness as a 
second site parameter. It defines the depth to a 
material (bedrock) with shear-wave velocity, vs 
of 1 km/sec or greater. 


E-ISSN: 2349 5359; P-ISSN: 2454-9967 
as the closest distance to the rupture plane 


2.1.2. Ground Motion Prediction Equations 


Ground motion prediction equations are used 
to measure ground-motion intensity, such as 
peak ground acceleration and _ spectral 
acceleration at the individual vibration period, 
including the characteristics of source, path and 
site parameters. Ground motion prediction 
equations are constructed by _ statistical 
regression analysis of many strong motion 
records. These ground motion models are 
developed for three distinct different tectonic 
zones: shallow crustal earthquakes, subduction 
zones, and stable continental regions. Chiou and 
Young ground motion model was developed 
from shallow crustal earthquakes [19]. The 
ground motion model is represented in 
Equations 2 and 3. 


In (Yee, ) = Cy + [CraFRvi + CrpF mi + C7 (Zrori — 4)](1 — AS)) + [C19 + C7a(Zrori — 4) JAS; + C2 (Mj — 6) 


C2 — C3 
+ 


n 


In(1 + efn(m—Mi)) + cqln[Raupi + cosh{c, max(M; — cy, 9)}] 


Cy2 
+ (Cy, —C in( R? i +c? )+ Cc ——— R ii 
( 4a re) RUPij RB y1 cosh|max(M; = cya) RUPij 


+coFywijtanh 
9a 


n(y;;) =i (Veet, ) +, -min (in (=), 0) + b> {ets(min(Vssop 1130)-360) _ e3(1130-360)} In ( 


1130 


ee 
ds (1 cosh[d¢. aE) 3 cosh[0.15 .max(0,Z1 9—15)] 


where,M Moment magnitude 
Reup The closest distance to the rupture 
plane (km) 

Rx Site coordinate (km) measured 
perpendicular to the fault strike 
from the surface projection of the 
up-dip edge of the fault rupture, 
with the down-dip direction being 


positive. 
Fuw Hanging-wall flag: 1 for Rx=0 and 0 
for Rx=0 
) Fault dip angle 
Fry Reverse faulting flag: 1 for 30° < 


dX <150° (combined reverse and 
reverse oblique), 0 otherwise; A is 
the rake angle 

Fy Normal faulting flag: 1 for -120° < 


(ses) i JR ypij + 2? Tori 


Rrupij+0.001 


(2) 


Yrefij _) 


4 


+ n; + Ej (3) 


Ryp Joyner-Boore distance to the 
rupture plane (km) 


X <-60° (excludes normal-oblique), 
0 otherwise. 

AS Aftershock flag: 1 if the event is an 
aftershock, 0 otherwise 


Vs30 Average shear wave velocity in the 
top 30 m (m/s) 

ZroR Depth to top of rupture (km) 

Z 1.0 Depth to shear wave velocity of 1.0 
km/s (m). 


Martin & Perez [20] used log normal 
distribution to model engineering data. Huyse et 
al., [21] demonstrated that 10-6 is the smallest 
annual _ probability exceedance for low 
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probability ground motions by using lognormal 
distribution. In this paper, lognormal 


The annual rate of exceedanceA(PGA > x) 
computed by the sum of the rates of PGA > x 
from each source to compute hazard curve is 
shown in Figure 2.The probabilistic seismic 
hazard analysis (PSHA) resulted in a horizontal 
PGA of 0.13g for design basis earthquake. The 
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distribution exceedance varying from 0.05 to 
0.8 g with intervals of 0.05has been generated. 


PGA value of the design basis earthquake is 
within the range of recommendations stated on 
the seismic hazard map of Ethiopia for return 
period of 475 years by [22] which is in the range 
between 0.1g and 0.16g as shown in Figure 3. 


Hazard curve 


0.0000 


Figure 2.Hazard Curve for the specific site 
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Figure 3: Seismic hazard map of Ethiopia based on the GSHAP data 
for a return period of 475 years [22] 


2.1.3 Deaggregation 


The hazard curve gives the combined effect of 
all magnitudes and distances on the probability 
of exceeding a given ground motion level. To 
provide insight into what events are the most 
significant contributions to the seismic hazard 
to determine the maximum _ considered 
earthquake (PGA), the hazard curve is broken 


down into its contributions from several 
sources for each related magnitude and source- 
to-site distance; which is the process of 
deaggregation. 


Deaggregation identifies the source, distance, 
and ground motion percentiles that contribute 
most strongly to the computed hazard. 
Deaggregation is obtained by dividing the total 
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hazard contributions based on distance and 
magnitude for the rate of an earthquake by the 
rate of all earthquakes [23,24]. The most 
common form of deaggregation is a two- 
dimensional (2D) deaggregation in magnitude 
and distance bins. In this method, the magnitude 
range is discretized within an equal space of 0.1. 
The deaggregation resulted in the most 
dominant magnitude of earthquake of 5.5 and 
5.7 and source to site distance of 20.23 km and 
15.27 km, respectively. This pair of values is 
considered as the most likely scenario which 
corresponds to the DBE for the specific site. 


2.1.4 Ground Motion Selection 


The selection and scaling of earthquake 
ground motions has a critical role in seismic 
load definition that will be applied to a structure 
during dynamic analysis. Since there are no 
accelerogram records for Ethiopian 
earthquakes, other records should be adapted, 
for which ground motion selection is required 
[25]. Ground motion selection — utilizes 
deaggregation results of magnitude and 
distance to identify causal events for a given 
target spectra value associated with an annual 
rate of exceedance [26]. 


To select acceleration time history from the 
database of Pacific Earthquake Engineering 
Research Center (PEER); the input factors are 
the pair of most dominant magnitude and 
distance from the deaggregation calculation, the 
top 30 m_ shear 


1.5 


Acceleration (m/s?) 
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velocity and fault mechanism of the specified 
area. Acceleration time history has been 
searched for an earthquake magnitude - 
distance combination range of 5.5-5.7 
magnitude and 15-30 km distance. Top 30 m 
shear wave velocity of 557 m/s has been used 
from the geotechnical report as an input on 
PEER ground motion database. The fault type 
was selected to be normal fault due to its 
dominance. 


Using the above PGA values, “Umbria Marche” 
earthquake has been selected from the database 
using the above mentioned inputs to 
characterize the area. The PGA value of this 
selected earthquake (0.163g) is used to scale the 
calculated PGA from PSHA value 0.13g. The 
scaling factor of the PGA becomes: 

0.13 


—— = 0.797 
0.163 


The most commonly employed ground motion 
scaling method involves multiplying all of the 
acceleration values of the time-acceleration 
pairs by a scalar value. This time domain 
scaling modifies the amplitude of the 
accelerations without affecting the frequency 
content or phasing. 


Z = 0.8 * 0.797 * 9.81 wh = 6.25 
& 


The scaled time acceleration histories are 
shown in Figure 4. 


Time (s) 


Figure 4: Scaled acceleration time history 


2.2 Three-dimensional finite element model 


The three dimensional Finite Element Method 
(FEM) software ABAQUS is used to model the 


pile- soil system under seismic excitation. The 
pile and the soil were modeled using linear 
brick elements with  eight-nodes and 
incompatible modes. A three dimensional 
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surface to surface contact element was used 
between the pile and the surrounding soil. Quiet 
boundaries were used to avoid wave reflection 
under seismic excitation by employing infinite 
elements. The quite boundaries do not transmit 
waves from near field to far-field without 
reflections. 


2.2.1 Soil Geometry and Properties 


The multilayered silty clay soil was modeled 
using Mohr-Coulomb constitutive model, with 
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the soil parameters obtained from _ the 
geotechnical investigation report of the project, 
summarized in Table 2. 


2.2.2 Pile Geometry and Properties 


The 15 m long concrete made pile is having 
circular cross sections of 0.6 m and 0.8 m 
diameters. The linear’ elastic material 
parameters of the pile are tabulated below 
(Table 1). 


Table 1: Geometry and properties of pile 


Unit Weight 


Young’s modulus 


Poisson’s ratio 


Pile diameter 


Pile length 


Pile head 


kN /m3 25 
kPa 30E6 
[-] 0.15 
m 0.6 and 0.8 
m 15 
[-] Free 


Table 2: Soil properties of different layers 


Symbol | t y E q’ c’ Ko 
Unit m kN/m3 | MPa ° kPa [-] 
Layer 1 | 4.5 17 8 0.3 18 12 0.691 
Layer 2 | 6 19 15 0.3 22 15 0.658 
Layer 3 | 4.5 20 35 0.3 24 20 0.626 
Layer 4 | 15 21 40 0.3 26 25 0.562 


Linear brick elements with eight-nodes of 
incompatible modes (C3D8I) were used to 
improve the computational time. The soil was 
assumed to be in dry state to avoid liquefaction 
problem. The dimensions of the soil shown in 
Figure Swere adopted from the 
recommendations of [27]. Material damping in 
soils was considered to be succeeded mainly 
through viscous damping. Therefore, 
traditionally, when computing material 
damping in soils, mass proportional damping is 
neglected and damping of the soil is achieved 
through _ stiffness proportional material 
damping. The governing equation of the system 
is given by: 


[M fu"}+[chu'}+[Kt+ {Fo} 


where: 


{u’}, {u’} and {u} are the acceleration, velocity, 
and displacement respectively. 


[M], [C] and [K] are mass, damping, and stiffness 
matrices, respectively. 


F(t) is force function. 
The damping matrix [C] = B[K], in which the 


damping coefficient 6 = oe , with, wa beeing the 
predominant circular frequency of loading, and 
t is the material damping ratio of soil assumed 
to be 5 % [10]. Predominant frequency is 
obtained from a Fourier spectrum drawn for the 
input wave from the analyzed time history as 
shown in Figure 6. 


2.2.3 Soil-Pile Interface 
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The contact at the circumferential pile 
boundary is modeled by defining the surface to 
surface contact between the pile and the 
surrounding soil. The pile-soil interaction is 
defined by tangential and normal contact 
behavior. Normal contact behavior used hard 
contact with penalty constraint enforcement 
method that allows any pressure to be 
transmitted between surfaces if they are in 
contact. Regarding the tangential contact, the 
shear behavior between two surfaces is defined 
by the Coulomb frictional law. The Coulomb 
friction model states that two contacting 
surfaces can carry shear stresses up to a certain 
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magnitude across their interface before they 
start sliding relative to one another [28]. The 
Coulomb friction model defines this critical 
shear stress, at which sliding of the surfaces 
starts as a fraction of the contact pressure, , 
between the surface, as given in Equation 5. 


C crit -_ Hp (5) 


Where tcrit the critical shear stress at contact 
surface, pis the coefficient of friction and, is the 
contact pressure between two surfaces. The 
friction coefficient was assumed to be isotropic 
and a value of p = 0.5 was used. 


Figure 5: Three dimensional view of the model 
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Figure 6: Fourier transformation of the earthquakes 


2.2.4 Numerical Accuracy and Stability 


The numerical simulation of dynamic Soil- 
Structure-Interaction (SSI) is controlled by two 
main parameters. These are the spacing of the 
nodes of the finite element model Ah and the 
length of the time step At. assuming that the 
numerical method converges toward the exact 
solution as At and Ah approach zero, the desired 
accuracy of the solution can be obtained as long as 
sufficient computational resources are available. 
In order to represent a traveling wave of a given 
frequency accurately, about 10 nodes per wave 


length are required. Fewer than 10 nodes can lead 
to numerical damping as the discretization misses 
certain peaks of the wave [29]. The value of the 
highest relevant frequency is 8.5 from Fourier 
amplitude analysis of input motion but 
recommended a frequency of 10 Hz to determine 
the appropriate maximum grid spacing for seismic 
analysis [29]. Thus, the maximum relevant 
frequency in this research is on acceptable range 
as shown Figure 6. 


The recommended formula relating grid spacing 
(Ah), time step of analysis (At), smallest wave 
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velocity of soil medium (v) and the highest 
relevant frequency (fma,) that is present in the 
model need to be found by performing a Fourier 
analysis of the input motion as given in equation 6 
and 7. 


-¥ 
ee (6) 
a= 
, (7) 


Where v is the highest wave velocity 


The smallest shear wave velocity of 110 m/s 
obtained from the geotechnical reports has been 
used in this research. 


2.2.5 Element Mesh Size 


For nonlinear numerical model of wave 
propagation problems, the time step At has to 
be limited for two reasons; i.e., accuracy and 
stability for all modes. If the time step in the 
finite element analysis is too large, the wave 
front can reach two consecutive elements at the 
same moment. This would violate a fundamental 
property of wave propagation and can lead to 
instability. The time step, therefore, needs to be 
limited. In this model, the smallest grid space Ah 
was 1.29 m and the maximum time step was 
0.012 s. 


Minimum mesh sizes that give accurate 
results and remove excessive calculation time 
and high computer memory requirements 
should be calculated. Accordingly, the mesh 
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sizes of the finite element calculation resulted in 
a size of 0.4 m. 


2.2.6 Boundary Conditions 


The important aspect of boundary conditions 
on dynamic analysis is to avoid wave reflection. 
In this paper, the soil was modeled into two 
regions, near field region consist finite element 
and far-field region consist infinite element. 
Infinite elements are used in boundary value 
problems with unbounded boundaries (infinite 
medium). During the dynamic analysis, infinite 
elements will provide “quiet” boundaries at the 
finite element model boundaries based on the 
model of [30]. The infinite element is not 
available in element definition on ABAQUS. 
Thus, the far-field region is first defined by 
three-dimensional acoustic elements with 8 
node linear brick elements (AC3D8). These 
AC3D8 elements are then replaced with CIN3D8 
elements, 8 noded three-dimensional linear 
brick elements. The advantage of using a 
combination of finite and infinite elements is 
minimizing computational efforts by reducing 
the volume of elements used in the analyses The 
infinite element boundaries also serve the 
purpose of boundary conditions, and no 
external boundary conditions need to be 
specified, as stated in the ABAQUS 
Documentation Manual. Regarding the 


combination of finite and infinite element near- 
field region and far-field region is shown in 
Figure 7a and the finite element discretization is 
shown in Figure 7b. 


Figure 7: a) Near and far field regions (Half model) _ b) Finite and infinite element discretization 


2.2.7 Loading Steps 


The computation of Finite Element Analysis is 
achieved by the following two loading steps; 
namely initial and dynamic loading steps. At the 
initial step of the model, actual in-situ 


conditions was simulated as an initial loading 
condition (in-situ stress conditions) before any 
dynamic or static external load. The geometry 
and properties of the soil, as well as the 
boundary conditions and the _ interaction 
between the soil and pile are defined in this 
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step, while the dynamic loading is modeled in a 
separate load step. 


In dynamic analysis, seismic load can be 
applied either as a displacement, velocity, or 
acceleration time history at the base of the 
model. The proposed design horizontal 
component of seismic acceleration time history 
corresponding to 475 years return period for 
Wolaita Sodo site is shown in Figure 4. 


Dynamic loading step is applied after geostatic 
load step. Seismic load is applied at the base of 
the soil-pile system in the form of acceleration 
time history after the selection and scaling of 
the designed time history. The seismic load is 
applied in both horizontal and_ vertical 
directions as the foundation supports a 
lightweight superstructure. The amplitude of 
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the vertical component is assumed to be 2/3 of 
the horizontal component (31). 


3. RESULTS AND DISCUSSION 


To determine the effect of pile size on pile 
head response, two pile diameters (0.6 m and 
0.8 m), which are practically used within the 
project premises have been considered while 
keeping pile properties, pile element, mesh 
sizes, and pile length constant. 


3.1 Pile Head Response 


Figure 8 shows the pile head response of 
displacement with dynamic time while Figure 9 
shows the acceleration response with dynamic 
time at pile head for 0.6 m and 0.8 m diameter 
pile, respectively, the pile head response is 
studied for both diameters of pile respectively. 


(b) 


= 
co] 
= 


= 
i—] 
id 
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Response (m) 
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nA 


-0.02 Time (s) 


| 
Input \) 


---- Pile Head Response 
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-0.08 


Figure 8: Lateral displacement under kinematic interaction a) diameter 0.6 m b) diameter 
0.8m 


It can be observed that the responses of the pile 
heads follow the pattern of the input acceleration. 
The input motions were amplified on the response 
diagrams, both for the acceleration and 
displacement. The maximum lateral deflections of 
the piles are 72.02 mm and 63.7 mm for 0.6 m and 
0.8 m diameter piles, respectively. Bouzid [32] 
stated allowable head deflection to be up to 5% of 
the pile diameter. Accordingly, the obtained pile 
deflections in this research exceed the acceptable 
range, the thresholds being 30 and 40 mm for 
diameter 0.6 m and 0.8 m respectively. This may 
be attributed to the ignored reinforcement that 
was not included in the computations to minimize 
computational efforts. 


It can be observed that the responses of the pile 
heads follow the pattern of the input acceleration. 
The input motions were amplified on the response 
diagrams, both for the acceleration and 
displacement. The maximum lateral deflections of 
the piles are 72.02 mm and 63.7 mm for 0.6 m and 
0.8 m diameter piles, respectively. Bouzid [32] 
stated allowable head deflection to be up to 5% of 
the pile diameter. Accordingly, the obtained pile 
deflections in this research exceed the acceptable 
range, the thresholds being 30 and 40 mm for 
diameter 0.6 m and 0.8 m respectively. This may 
be attributed to the ignored reinforcement that 
was not included in the computations to minimize 
computational efforts. 
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Figure 9: Acceleration response at pile head of diameter 0.6 m and 0.8 m piles, respectively, with 
dynamic time under kinematic interaction 


3.2 Embedment Depth Effect 


The lateral deformation of the pile at different 
depths for both pile diameters under 
consideration are shown in Figure 10. The 
lateral displacement of the pile decreases with 
increasing embedment depth, due to the 
increasing soil stiffness with dept. The bending 
moment diagram along the pile length is shown 
in Figure 11 for both pile diameters. The 


maximum bending moment is attained at the 
bottom of the pile. The pile diameter affects the 
amplitude of the bending moments close to the 
bottom of pile, the higher moment being 
exhibited for the larger pile diameter, due to the 
increased stiffness. The patterns of bending 
moment under dynamic load for both diameters 
are the same and this result is similar to 
previous literature [33]. 
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Figure 10. Lateral displacement under kinematic interaction a) for diameter 0.6 m pile 
b) for diameter 0.8 m pile 
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Figure 11: Bending moment along the pile length a) Diameter 0.6 m b) Diameter 0.8 m 


4. CONCLUSIONS 


During ground shaking, pile foundations are 
subjected to lateral forces. Even though the 
primary function of a pile foundation is to carry 
vertical loads from the superstructure, it has to 
withstand lateral forces due to seismic actions. 
The seismic response of the soil is affected by an 
interaction between the structure’ and 
surrounding soil. 


In this research, seismic hazard analysis for 
the specified electric converter station area was 
carried out. Probabilistic seismic hazard 
analysis has been carried out for design basis 
earthquake, corresponding to 475 years of 
return period and its peak ground acceleration 
is determined to be 0.13¢ that falls in the range 
of PGA values determined by previous 
researchers. 


The response of the pile under seismic load in 
layered soil has been carried out by employing 
finite element analysis software ABAQUS. The 
dynamic analyses results by considering plain 
concrete sections of the pile showed that the 
displacement of the pile is not within 
permissible ranges of the literature, which calls 
for further detailed analyses. Otherwise the use 
of the smaller diameter pile shall be 
disregarded, as the calculated displacement is 
more than twice the allowable. Qualitatively, the 
results of current investigation are in good 
agreement with the findings of previous 
researches and theoretical expectations. 
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